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Using Random Phase Approximation, we show that a crossed-chains model of a spin- 1/2 Heisenberg 
spins, with frustrated interchain couplings, has a non-dimerized spin- liquid ground state in 2D, with 
deconfined spinons as the elementary excitations. The results are confirmed by a bosonization 
study, which shows that the system is an example of a 'sliding Luttinger liquid'. In an external 
field, the system develops an incommensurate field-induced long range order with a finite transition 
temperature. 
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Geometrically frustrated magnets have attracted 
strong interest in recent years, primarily because they 
provide us with a potentially direct route to the ac- 
tively sought "spin liquid" phases of two- (2d) and three- 
dimensional (3d) strongly correlated electron systems. 
Begun in the 50's this line of research was dramat- 
ically revitalized by Anderson's proposal ||^ that such 
a liquid of spin singlets is, upon hole doping, adiabat- 
ically connected to the ground state of layered cuprate 
superconductors. Whether or not this is so remains to 
be seen, but the following "question of principle" — can 
the ground state of a 2d or 3d magnetic system be liquid 
like? — has become one of the most debated questions 
in condensed matter physics. Significant insight into this 
problem has been gained recently via the dual description 
of frustrated antiferromagnets in terms of gauge Ising 
models [^^. In this paper we, however, take a more 
direct, experimentally- motivated approach, following re- 
cent work by Bocquet et al. and an earlier paper by 
two of us and P. J. Freitas ||^. 

The "spin liquid" (also known as "resonating valence 
bond", or RVB, phase) is defined as a liquid of singlet 
spin pairs covering the lattice, and is characterized by 
the absence of long-range order (LRO), unbroken spin 
rotational and translational symmetries, and elementary 
excitations with fractional spin 1/2 {spinous). Its phys- 
ical relevance has been highlighted by the experimental 
observation of spin liquid-like behavior in the kagome 
lattice compound SrCr8_2;Ga4+2;Oi9 j^] and several py- 
rochlore antiferromagnets |l^, e.g. CsNiCrFg. Both ma- 
terials are based on frustrated units — triangles (kagome) 
and tetrahedra (pyrochlore) — combined in site shar- 
ing arrangements, which leave the spins in the classical 
ground-states highly underconstrained pd| . 

On the other hand, it is well known that a spin liq- 
uid state is realized in the antiferromagnetic Heisenberg 
chain (HAFC). It is then natural to try to build up 2d 
(or 3d) spin liquid state from HAFCs - the theoretical in- 
vestigation of this possibility is the subject of this paper. 

Elementary excitations of the single chain are decon- 
fined gapless spin- 1/2 spinons, which can be visualized 



as domain walls separating domains of different orienta- 
tion (up and down) of staggered magnetization n{x) = 
(— l)^S'(x) along the chain. Consider now a parallel array 
of HAFCs coupled by unfrustrated exchange J± in the di- 
rection transverse to the chain: no matter how small J± 
is, it immediately leads to the confinement of spinons be- 
cause the energy of two domain walls grows linearly with 
separation between them. Fig. la. This is an intuitive 
reason, confirmed by detailed calculations for the 
stabilization of LRO and appearance of spin-1 magnons 
(which are bound states of two spinons) in this situa- 
tion. It is clear, however, that this argument fails if each 
spin is coupled, in a transverse direction, to the zero- 
spin combination of spins on neighboring chains, i.e. if 
the transverse coupling is frustrating. Motivated by this 
simple argument and recent experiments on CS2CUCI4 
[p^ and, especially, Na2Ti2Sb20 |lj], we investigate the 
crossed-chains model Q (CCM) — the two-dimensional 
Heisenberg model on the lattice shown in Fig. lb, 



Hcc — J ^ ^ SiSj -l- J ^ ^ SiSj 



(1) 




The intra-chain exchange J couples neighboring spins 



FIG. 1. (a) Unfrustrated array of HAFCs. Arrows repre- 
sent staggered magnetization n{x) . Solid (dashed) lines stand 
for exchange J (J±), whereas strings mark broken J_l bonds. 
Domain walls are indicated by crosses, (b) The crossed-chains 
lattice. Filled circles (squares) indicate spins on the vertical 
(horizontal) chains. Solid (dashed) lines denote intra-chain J 
(inter-chain J') exchange interactions. 

in a given row (column) which form a square mesh of 
HAFC along the X and Y axes, whereas the inter-chain 
J' is the nearest neighbor coupling between spins in rows 
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and columns. The crossed-chain Hamiltonian interpo- 
lates between decoupled HAFCs (J' = 0), the 2d py- 
rochlore lattice [J' = J), and 45°-rotated square lattice 
(J — 0) Notice that, when neighboring spins are an- 
tiparallel along the chains, the inter-chain contribution to 
the effective field acting on a given spin is zero. 

1. RPA analysis of the crossed- chains model in the 
weak-coupling limit J' <C J. The defining feature of 
the HAFC is the absence of coherent spin-1 magnon ex- 
citations, already mentioned above. The spectrum of 
magnetic excitations consists of spin- 1/2 noninteracting 
spinous, and the usual spin-1 magnon is an incoherent 
two-spinon excitation. It is then natural to use an ap- 
proach which properly accounts for this important one- 
dimensional feature in the "decoupled" limits J' / J 
and/or T ^ J' . This is provided by a Random-Phase- 
Approximation (RPA) 17 which has the meaning 



of an expansion in the inverse coordination number of 
the lattice z± {— 4 for CCM). Applied to the crossed- 
chains lattice of Fig. lb, it gives the following result for 
the dynamical spin susceptibility of CCM 

XB^PA{!^,K-,ky) = {xi{^,K) +xi{^,ky) - 2J'{k) x 

Xi{oj,kx)xi{^.ky)){l-{J\k)fxi{^.kx)xi{uJ.ky)) ' (2) 



— -I- 7rtanh(27ra;o) - 2Im*(i -I- ixo) = 0. (5) 
xq 4 

This equation is very similar to the one obtained for the 
spin model describing CS2CUCI4 - in that case the co- 
efficient of the first term in (|5|) is 2 times smaller. This 
minor difference is, however, extremely important: unlike 
the CS2CUCI4 model, Eq.(^ has no solution. Thus, RPA 
analysis predicts no ordering instability down to, and in- 
cluding, T = 0: the crossed chains remain decoupled and 
provide us with an example of two-dimensional spin liq- 
uid with deconfined spinous as elementary excitations. 
The ineffectiveness of J' in destabilizing Id behavior of 
the crossed-chains is clearly seen from Fig.^ where we 
compare xrpa (^, fc. A:) (2 ) with twice the susceptibility of 
a single chain xi(a', fc) (3). Even for unrealistically large 
J' = J the difference between these two susceptibilities 
is hardly observable. The same is true for corresponding 
structure factors Sy{uj, k) = —(l—e~^/'^)~^lviiXv{'^, k, k) 
[v — RPA, 1) which arc compared in Fig.|^. 



where momentum k — {kx,ky) is measured from the an- 
tiferromagnetic momentum (tt, tt) (lattice spacing is set 
to unity). The dynamical susceptibility of a single hor- 
izontal {kx) and vertical {ky) chain, xi(tt',A:), is known 
exactly (for details see, e.g., 0), 



Xii^ik) 



vk , 



Pri(x) = 



2(27r)3/2T 
^-ix) 



T{1-1- ix) 



(3) 



where r(x) is the Gamma function, A = 2A.27J is 
the high-energy cutoff, numerical value of the pre-factor 
was calculated in ||l^, rj = 1 for the isotropic HAFC, and 
V = TT J/2 is the spinon velocity. 

The specific lattice structure of CCM is encoded in the 
Fourier transform of the inter-chain exchange interaction 



J'(fc) = 2J'sin(fc^/2)sin(fcj^/2), 



(4) 



and its frustrating character is clear from the fact that 
J'(0,0) = 0. 

An ordering instability, if any, should show up as a 
divergence in Xi?.PA (0,fco) at some critical temperature 
To and momentum fco = (fco,±fco) along two diagonal 
directions where antiferromagnetic fluctuations are the 
strongest. Extremizing the denominator of (j^) with re- 
spect to fco one arrives at the following implicit equation 
for the ratio xq = fco/(47rro) {^(x) is the digamma func- 
tion) 




FIG. 2. —\mxRPA(oJ, k, k) (solid line) and —2 x Imxi(ii', k) 
(dashed line) vs. k for J' = J, lj = O.IJ, T = O.OIJ. 




FIG. 3. SRPA{Lo,k) (solid line) and 2 x 51(1.^,^) (dashed 
line) vs. u>/J. k — 0.06 and other parameters as in Fig.^. 

The reason for this is purely geometrical and is best 
illustrated at T = 0: the 1/k divergence of Xii'^ = 0,k), 
Eq.(||), is suppressed by the fc^ factor from J' (A;), Eq.(||). 
This, again, should be contrasted with the CS2CUCI4 
lattice where the inter-chain exchange scales as k and 
fails to suppress a remaining weakly diverging factor 
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v/ln(l/|fc|) due to a marginally irrelevant Umklapp in- 
teraction, which eventually does lead to the instability at 
some incommensurate wavevector . 

Going back to the crossed-chains model it should be 
noted that the ordering instability may, in principle, hap- 
pen at some "noncritical" fco in the bulk of the Bril- 
louin zone (BZ) when the ratio J' / J exceeds some criti- 
cal value. Analysis of this interesting possibility requires 
knowledge of xi(0, fc) for all values of momentum inside 
the BZ, which is currently not available. However, there 
exist a closely related model for which the static suscep- 
tibility is known exactly for all momenta - the Haldane- 
Shastry chain (HSC) ||^. At T = its susceptibility 
diverges as 1/fc near the antiferromagnetic momentum 
TT, similar to the HAFC, and remains finite everywhere 
else, XHs{^,k) ~ arccos |fc|/(27rw|fc|v'l — k"^)- Replacing 
Xi(0, k) in Eq.(||) with x_f/s(0, k) in order to get an idea 
of what may happen in the CCM as the ratio J'/ J in- 
creases, we find an instability at fco ~ 0.6477r when the 
ratio of exchanges exceeds (J'/J)crit — (7i'/2)'^ ~ 3.87. 
This high value of J^.^.-^^ in the HS model on the crossed- 
chains lattice gives additional support to the robustness 
of the spin liquid state. 

2. Effect of magnetic field. Interesting behavior is ex- 
pected to develop when the crossed-chains model is sub- 
jected to a uniform magnetic field H . An applied field 
breaks spin-rotational symmetry between transverse and 
longitudinal susceptibilities 0. The transverse one re- 
mains dominated by excitations from the gapless states 
at fc = whereas the maximum of the longitudinal 
susceptibility Xi^ shifts to the incommensurate points 
k± = ±6{H), SiH) = 2ttM{H), where M{H) is the 
magnetization of a single HAFC. Since J'(fc) is not af- 
fected by the field, the divergence of Xi^(0, fc k± ) is 
not compensated by inter-chain exchange anymore and 
an ordering instability does develop. It signals a transi- 
tion into an incommensurate field-induced ordered phase 
(FILRO). Both the ordering temperature Tq and the or- 
dering momentum fco can be calculated from a properly 
generahzed Eq.(||). Repeating steps that led to (||) we 
obtain (fco is measured from k± now) 



2^ro ..fcoi'^, ^ 

— cot(— — ) + 



TT sinh(27r2;o) 



cosh(27ra;o) — cos{^) 



-2Im*(- +i2;o) = 0. 
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(6) 



This equation has to be solved simultaneously with an 
equation for To, which is simply the condition that de- 
nominator of XRPAi^jko) is equal to zero. Notice that 
both the scaling exponent ri{H) and the spinon ve- 
locity v{H) decrease with magnetic field, in particular 
ri{H) = 1/2 for saturating field 7J > 2J j?). How- 
ever, analytical solution is still possible in the limit 
of weak applied field when 5{H) thanks to the 
following inequalities vkg ^ Tq, fco <C (5 which can 



be checked a posteriori. We find incommensurate or- 
dering with fco,± = ±ciJ(-^''+-^^/(^''~^^ developing at 
To = C2v{Sy^J'/vf'^^'^^'^-^\ where Ci,2 are weakly in- 
dependent constants of order 1. So that for if ^ the 
scaling is fco,± ~ ±S^, Tq ~ J'6'^. Notice that 2d ordering 
occurs at momenta different from k± at which suscepti- 
bility has a maximum for independent chains, although 
the difference is probably too small to be observed ex- 
perimentally. Observation of such a FILRO in materials 
such as Na2Ti2Sb20 would be a direct signature for the 
applicability of the coupled-chains model presented here. 

3. Scaling analysis: sliding Luttinger liquid phase. 

We now employ bosonization to investigate effect of the 
inter-chain interaction V at finite temperature T . We 
consider the mesh made of equal number of horizontal 
{h) and vertical {v) chains, Nh = Ny = N , each of length 
L = Na, where a is the lattice spacing. In the continuum 
limit each site spin is represented by a sum of uniform J 
and staggered n fields Q. In addition, J^^^ = J^y^ ^ -|- 
^u/Zi L' where J'^/j^ are chiral WZW currents. The inter- 
chain interaction reads 



V ^ g dxdyJhix, y) ■ Jy{x, y) 



(7) 



where g = 23' . Note that due to the geometry of the 
problem the staggered magnetization n, which has scal- 
ing dimension A(r?) = 1/2, does not show up in this 
equation. As a result, the scaling dimension of the in- 
tegrand in is 2. The bare Hamiltonian describes 
total of 27V independent vertical and horizontal chains. 
From here follows the defining feature of our model: J —J 
correlations are non-zero only for like currents (i.e. h-h 
and v-v) of the same chirality (i.e. R-R and L-L) which 
belong to the same chain, e.g. 



(J„%(x,y,r)j;^ (0,0,0)) = aS{x)S^,,,i7rT/v)^ 



''87r2^sin2 



ttT{vt + iy)/v] 



sm' 



[ttT(vt — iy)/v 



} (8) 



and (JhJh) is obtained by replacing u — s- ft,, y x. 
The 'same-chain' condition is contained in an important 
delta-function which appears with pre-factor a because 
5{x — x') — 5nn'/a. Conformal invariance of was 
used to write (|) at finite T HJI^. 

Our idea is to consider corrections to the bare free en- 
ergy of 2N chains Fo/L'^ = TTv/{2a^) - irT'^/iiva) (mo- 
mentum cut-off A ~ n/a was used to get the first term). 
First correction is — —T(V^)/2 where 



(V^) = igaf / dxdy 



l/T 



drdr'iJhix, y, T)Jh{x, y, r')) 



x{Jy{x,y,T)Jy{x,y,T')) 



(9) 



Notice that due to delta-functions in (|^) spacial coordi- 
nates of currents in (|^) are forced to coincide. Short- 
distance divergence of the integrand in (j9|) is regularized 
by a, leading to 
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c^) - (^) 



2tt 



d(j) 



27r[cos0 - cosh(27rra/u)]2 



(10) 



The resultant correction to the free energy is irrelevant, 



' 2a3 V87r2i; 



\\-U:i^Talvf}. (11) 



This can be understood as follows. In the weak-coupling 
limit, where each crossing can be treated independently 
from others, interaction V acts at a 'point rather than 
along the entire length of two chains forming that partic- 
ular crossing. This is the physical reason for the absence 
of integration over the relative spatial coordinates in (|^). 
It is well known that for the point interaction marginal 
dimension separating relevant from irrelevant perturba- 
tions is 1 |23j. Since A(V^) = 2 > 1, F is irrelevant. 

At this order the ground state of the CCM is that of 
decoupled spin chains. At T = there are no correla- 
tions between spins on different chains due to irrelevancy 
of V but at finite T weak inter-chain correlations with 
a square C^y symmetry will be present [ p4[ , in qualita- 
tive agreement with our RPA expression (^. Dimen- 
sional estimate of higher order in g corrections show that 
their contribution is at most marginal. This allows us to 
identify the CCM as a SU (2) invariant spin liquid with 
deconfined spinous as elementary excitations. Following 
earlier works on models with U{1) symmetry [^|4| it 
can be called a sliding Luttinger liquid. 
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Note added: Since the submission of our work nu- 
merical studies have further addressed the J' = J limit 
of the model (the 2d pyrochlore lattice) |2^ . They find 
a valence-bond crystal phase with a spin-gap. These re- 
sults are not in contradiction with our findings. In fact, 
making the reasonable assumption that the ground state 
energy in the valence-bond crystal phase and the sliding 
Luttinger liquid phase is weakly dependent on J'/ J, leads 
to an estimate for a quantum phase transition between 
the two at J' / J = 0.85. This allows for a wide range 
< J'/ J < 0.85 of values where the sliding Luttinger 
liquid phase could exist. 
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